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ABSTRACT 

We have derived new abundances of the rare-earth elements Pr, Dy, Tm, Yb, 
and Lu for the solar photosphere and for five very metal-poor, neutron- capture 
r-pro cess- rich giant stars. The photospheric values for all five elements are in 
good agreement with meteoritic abundances. For the low metallicity sample, 
these abundances have been combined with new Ce abundances from a com- 
panion paper, and reconsideration of a few other elements in individual stars, 
to produce internally-consistent Ba, rare-earth, and Hf (56 < Z < 72) element 
distributions. These have been used in a critical comparison between stellar and 
solar r-process abundance mixes. 

Subject headings: atomic data - Sun: abundances - stars: abundances - stars: 
Population II - stars: individual (CS 22829-052, CS 31082-001, HD 115444, HD 
221170, BD+17 3248) 



1. INTRODUCTION 

Early Galactic nucleosynthesis studies have been invigorated over the last decade by 
the discovery of many low metallicity halo stars with abundance distributions that depart 
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significantly from that of our Solar System. The neutron-capture elements (Z > 30, hereafter 
n-capture) as a group exhibit particularly large star-to-star abundance variations with re- 
spect to Fe-peak elements. For example, data from a number of surveys collected in Sneden, 
Cowan, Sz Gallino (2008) show an abundance range in the rare-earth element Eu of at least 
-0.5 < [Eu/Fe] < +2.0 at metallicities [Fe/H] < -2.50; see their Figure 14. 

The n-capture abundances in the solar system and in most metal-rich Galactic disk 
stars arise from the combined effects of prior rapid and slow n-capture synthesis events 
(the "r-process" and "s-process", respectively). The n-capture abundance patterns in low 
metallicity stars, however, vary widely. Examples have been found with element distributions 
that are consistent with the r-process, the s-process, and a variety of mixes in between 
these two extremes. These stars are thus natural test cases for n-capture nucleosynthesis 
predictions. 

Rigorous tests of r-process and s-process theories require very accurate n-capture abun- 
dances in metal-poor stars. Good abundance determinations result from effort on all fronts: 
acquisition of very high resolution, low noise spectra of the stars; construction of realistic 
model stellar atmospheres; analysis of the spectra with few limiting simplifications; and 
improvement in basic atomic and molecular data. We have taken up the last considera- 
tion in the present series of papers: Lawler, Bonvallet, & Sneden (2001a), Lawler et al. 
(2001b), Lawler et al. (2001c), Den Hartog et al. (2003), Lawler, Sneden, & Cowan (2004), 
Den Hartog et al. (2005), Lawler et al. (2006), Den Hartog et al. (2006), Lawler et al. (2007), 
Sobeck, Lawler, & Sneden (2007), Lawler et al. (2008b), and Lawler et al. (2009). We have 
concentrated most of our efforts on: (a) improving the basic laboratory data for (mostly) 
rare-earth ionized species that are detectable in metal-poor stars; (b) applying these data 
to derive new solar spectroscopic abundances and comparing these photospheric values to 
solar-system meteoritic data (Lodders 2003); and (c) extending the abundance analyses to 
a few well-studied low metallicity giants that are enriched in the products of the r-process. 
Our most recent study (Lawler et al. 2009) reports improved transition probabilities for 921 
lines of Ce II. The present paper culminates this series with new solar and stellar analyses of 
Pr, Dy, Tm, Yb, and Lu. These elements all have good laboratory studies of their first ions 
in the literature, but have not been systematically subjected to solar /stellar analyses in the 
same manner as have other rare earths. 

In this paper we expand the standard definition of the rare-earth elements from the 



1 We adopt the standard spectroscopic notation (Heifer, Wallerstein, & Greenstein 1959) that for elements 
A and B, [A/B] = logi (N A /N B )^ - logio(N A /N B )0. We use the definition log e(A) = logi (N A /N H ) + 
1 2.0, and equate metallicity with the stellar [Fe/H] value. 
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lanthanides (57 < Z < 71) to include two adjacent elements Ba (Z = 56) and Hf (Z = 72), 
and adopt the collective shorthand notation "RE" for them. This broad definition covers a 
contiguous set of elements that have similar properties for stellar spectroscopy. In particular, 
these elements have relatively low first ionization potentials, 5.2 eV < IP < 6.8 eV, and 
thus are almost completely ionized in the solar photosphere and in the atmospheres of low- 
metallicity giant stars. Their only detectable spectral features arise from their first ionized 
species. Element groups in the Periodic Table immediately preceding the REs (e.g., I, Xe, 
Cs) and following them (e.g., Ta, W, Re) have very different atomic properties. For various 
reasons traceable to very low abundances, Saha/Boltzmann energy level population effects, 
and/or lack of accessible transition wavelengths, these elements just outside the RE group 
are inaccessible to most stellar spectroscopic detection efforts. 

In §5] we review the solar and stellar spectroscopic data and outline the abundance 
derivation methods. Results for individual elements are given in §3j We summarize the total 
RE abundance sets for the solar photosphere and r-pro cess-rich metal-poor giant stars in 
£JD Finally, in S|5] we use the stellar RE abundance distributions in a critical examination of 
r-process predictions. 



2. SPECTROSCOPIC OBSERVATIONS, REDUCTIONS AND ANALYSES 

For most of our stars, we analyzed the same high-resolution spectra that have been 
used in previous papers of this series. Additional descriptions of these stellar spectra can 
be found in their original studies: BD+17 3248, Cowan et al. (2002); CS 22892-052, Sneden 
et al. (2003); HD 115444, Westin et al. (2000); HD 221170, Ivans et al. (2006); see also 
Cowan et al. (2005). The spectroscopic data sets employed in our analysis are summarized 
in Table [TJ For each of the instrumental setups listed, we report the useful wavelength range, 
and wavelength-dependent values of the signal-to- noise S/N, resolving power R, and quality 
factor per resolution element F (sometimes also referred to as figure-of-merit), at selected 
wavelengths X app . 

Data reduction for the Keck and McDonald data have been detailed in previous papers 
of this series and have largely relied on IRAF@ and FIGARO i For the recently acquired 



2 IRAF is distributed by NOAO, which is operated by AURA, under cooperative agreement with the 
NSF. 

3 FIGARO is a part of the "Starlink Project" , which is is now maintained and being further developed 
by the Joint Astronomy Centre, Hawaii. 
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Magellan/MIKE data, we employed the MIKE Pipeline software^] (Kelson et al. 2000, Kel- 
son 2003). All of the data received final pr ocessing including continuu m normalization 



and telluric feature removal using SPECTRE (IFitzpatrick fc Snedenl 119871 ). Finally, for the 



solar analyses we employed the very high resolution (R = > 300,000), very high signal-to- 
noise (S/N > 1000) center-of-disk photospheric spectrum of Delbouille, Roland, and Neven 
(1973 j| 

The abundance analyses used the same methods that have been described at length in 
previous papers of this series. Here we summarize the main points; the reader should consult 
Lawler et al. (2009) and references therein for details. 

For each species, we begin with computations of relative strengths of all lines, in order 
to trim the often extensive laboratory line lists to a set that might produce detectable 
absorption in the solar photosphere and in our program stars. Line absorption coefficients 
are proportional to products of oscillator strengths and absorber number densities. In a 
standard LTE analysis, Boltzmann/Saha statistics describe the populations of atoms in 
various ionization stages and electronic levels. As discussed in SJTJ the REs have low ionization 
potentials, and thus exist almost completely as singly ionized species. Saha corrections 
for other ionization states can be neglected. Therefore the relative strength factors for 
ionized-species RE elements can be approximated by log (eg/) — 0%, where e is the elemental 
abundance, gf is the transition probability, 9 = 5040/T is the reciprocal temperature, and 
X is the excitation energy. 

Almost all easily detectable RE lines are of low excitation, x ^ 1 e V, so the relative 
line strengths are not very sensitive to temperature. Choosing 9 = 1.0 as a rough mean 
of the solar and stellar reciprocal temperatures, and adopting approximate solar abundance 
values for each element under consideration, we computed relative strength factors for Pr II, 
Dy II, Tm II, and Lu II lines, using the laboratory data that will be discussed in the ap- 
propriate subsections of $3j We did not perform such computations for Yb II, as it has 
only two very strong lines of interest for abundance analyses (see §3.4p . The results of this 
exercise are displayed in Figure [H With horizontal lines we mark the approximate mini- 
mum relative strength value for lines that can be considered "strong". Such lines are those 
with evident saturation in their equivalent widths (EWs), which for the Sun empirically is 
\og(EW/X) ~ —5.3. We similarly mark the approximate strength value at which photo- 
spheric lines have \og(EW/X) ~ —6.5, too weak to be routinely used in abundance solar 



The MIKE Pipeline is available from the Carnegie Observatories Software Repository at 



http: / /www. ociw.edu/Code/mikc/ 



5 Available at http://bass2000.obspm.fr/solar_spect.php 
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abundance analyses. 

Figured] can be compared to similar plots for other RE elements in some of the previous 
papers of this series. Some general remarks apply to all RE ions. Most REs have complex 
energy structures, leading to large numbers of transitions. Their relative strength factors 
increase with decreasing wavelength; these usually are transitions from the lowest energy 
levels with the largest log gf values. The most fertile regime for RE transitions is the 
near-UV domain, A < 4000 A. Unfortunately, the strong-line density of all species increases 
in this wavelength range, and many promising RE transitions are hopelessly blended with 
(usually) Fe-peak lines. Finally, as is evident in Figure [U very few RE ions have detectable 
transitions in the yellow-red (A > 5000 A) spectral region of the solar spectrum. Comments 
on the line strengths of individual species will be given in §3J These same strength factors 
turn out to work reasonably well for the r-process-rich giant stars. Their combination of 
cooler temperatures, more extended atmospheres, metal poverty, and enhanced n-capture 
abundances yields line strengths that are similar to or somewhat larger than those for the 
Sun. 

We eliminated lines with relative strength factors that fell below the probable detection 
limits, and searched solar and stellar spectra for the remaining lines. In this effort we 
employed the large Kurucz (1998)0 line list, the solar line identifications of Moore, Minnaert, 
& Houtgast (1966), and the observed spectra described above. With these resources we were 
able to discard many additional lines that proved to be too weak and/or too blended to be 
of use either for the Sun or for the r-process-rich stars. 

We then constructed synthetic spectrum lists for small spectral regions (4-6 A) sur- 
rounding each promising candidate line. These lists were built beginning with the Kurucz 
(1998) atomic line database. We updated the n-capture species transition probabilities with 
results from this series of papers, including the laboratory data cited below for Pr, Dy, Tm, 
Yb, and Lu. We also used recently published log gf values for Cr I (Sobeck, Lawler, & 
Sneden 2007) and Zr II (Malcheva et al. 2006). Lines missing from the Kurucz database but 
listed in the laboratory studies or in the Moore et al. (1966) solar line atlas were added in. 
In spectral regions where molecular absorption is important, we used the Kurucz data for 
OH, NH, MgH, and CN, and Plez (private communication) data for CH. 

We iterated the transition probabilities through repeated trial spectrum syntheses of 
the solar photosphere (and sometimes one of the r-process-rich giant stars). For the Sun, as 
in previous papers of this series we adopted the Holweger & Muller (1974) empirical model 
photosphere, and computed the synthetic spectra with the current version of Sneden's (1973) 



6 Available at http://cfaku5.cfa.harvard.edu/ 
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LTE, 1-dimensional (ID) line analysis code MOOG. In these trial syntheses, no alterations 
were made to the lines with good laboratory logp/'s. On occasion, obvious absorptions 
without plausible lab or solar identifications were arbitrarily defined to be Fe I lines with 
excitation energies % = 3-5 eV and log gf values to match the photospheric absorption. 
We discarded all candidate RE lines that proved to be seriously blended with unidentified 
contaminants. 

Final solar abundances for each line were determined through matches between the 
Delbouille et al. (1973) photospheric center-of-disk spectra and the empirically-smoothed 
synthetic spectra. The same procedures were applied to the observed stellar spectra (Tabled]) 
and synthetic spectra generated with the model atmospheres whose parameters and their 
sources are given in Table [2J 

3. ABUNDANCES OF PR, DY, TM, YB, AND LU 

In this section we discuss our abundance determinations of elements Pr, Dy, Tm, Yb, and 
Lu in the Sun and the r-pro cess- rich stars. Tables |3llH and [5] contain the mean abundances 
in the solar photosphere and in the r-process-rich low metallicity giants for these elements 
and for other REs that have been analyzed in previous papers of this series. The full suite 
of elements will be discussed in §U Table [3] also gives estimates of r-process abundance 
components in solar-system meteoritic material. These data will be discussed in more detail 
in §3 

3.1. Praseodymium 

Pr {Z = 59) has one naturally-occurring isotope, 141 Pr. The Pr II spectrum has been 
well studied in the laboratory, with transition probabilities reported by Ivarsson et al. (2001; 
hereafter IvOl), Biemont et al. (2003) and Li et al. (2007; hereafter Li07), as well as numerous 
publications on its wide hyperfine structure. We will consider the hfs data in more detail in 
the Appendix. 

We adopted Li07 as our primary transition probability source. This is the most recent 
and largest set, 260 lines, of purely experimental measurements (Li07 combined their own 
branching fractions with previously published lifetimes). IvOl also conducted a smaller Pr II 
lab study, reporting log g /values for 31 lines. However, their list includes four lines not 
published by Li07. Therefore we considered both Li07 and IvOl data sets in our abundance 
determinations. In Figure [2] we plot the differences between individual IvOl and Li07 log gf 
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values, using different symbols to distinguish those lines employed in our abundance analyses 
from those that proved to be unsuitably weak or blended. There is generally good agreement: 
ignoring the five obviously discrepant lines that are labeled by wavelength in the figure, the 
mean difference is < log gf Iv01 — log gf Lm > = +0.03 ± 0.01 (a = 0.06, 23 lines). Comments 
on individual lines in common are given below. Biemont et al. (2003) also published logg/'s 
for 150 Pr II lines. However, their values were determined by combining experimental Pr II 
lifetimes and theoretical branching fractions, which are very difficult to compute for the 
complex RE atomic structures (e.g., Lawler et al. (2008a). 

Moore et al. (1966) give 21 Pr II identifications for the solar spectrum. However, most 
of them are very weak and/or blended. An early study by Biemont et al. (1979) has a 
good discussion of the benefits and disadvantages of many of these lines for photospheric 
abundance work. They used nine lines to derive log e(Pr) Q = 0.71 ± 0.080 with individual 
lines contributing to the average with different weights. Only three of these lines were 
considered to be high-quality ones. More recently, Ivarsson, Wahlgren, & Ludwig (2003) 
employed synthetic/observed spectral matches to suggest log e(Pr) = 0.4 ± 0.1, more than 
a factor of two smaller than the meteoritic value of log e(Pr) me t = 0.78 ± 0.03. 

We searched for useful Pr II lines in the solar spectrum by first identifying them in 
CS 31082-001, which is the most extreme r-process-rich metal-poor star of our sample: 
[Fe/H] = -2.9, [Eu/Fe] = +1.6 (Hill et al. 2002). This star's low metallicity and large 
[n-capture/Fe-peak] abundance ratios combine to yield many strong (and often essentially 
unblended) candidate transitions. Inspection of the CS 31082-001 spectrum yielded 43 lines 
from Li07 and an additional 3 lines from IvOl that merited abundance consideration. Prelim- 
inary synthetic spectrum calculations suggested that 13 of these candidate lines were either 
too weak or too blended in both CS 31082-001 and the Sun. The wide hyperfine structure 
of all prominent Pr II lines made this exercise much easier than it would be in searches for 
lines with no hfs. In Figure [3] we illustrate this point with synthetic and observed spectra of 
the strong 4408.8 and 4179.4 A transitions. Visual inspection of the Pr II profiles suggests 
that their full-width-half-maxima are FWHM ~ 0.4 A, while observed and synthetic profiles 
of single lines (e.g., 4178.86 A Fe II and 4179.59 A Nd II) have FWHM ~ 0.25 A 

Wavelengths of the remaining useful Pr II lines are given in Table El along with their 
excitation energies and the Li07 and IvOl transition probabilities. In Figure [2] one sees 
five lines with large log gf discrepancies between these studies. Three of the lines were not 
involved for our abundance studies and so we cannot comment further on them. Li07 caution 



7 Throughout this paper we will use the subscript symbol to indicate solar photospheric values, and 
the subscript met to indicate solar system meteoritic values from Lodders (2003). 
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that 5219.1 A is blended on their spectra. We adopted the IvOl value for this line. Finally, 
the difference between IvOl and Li07 for 5322.8 A is 0.2 dex, but abundances derived with 
the Li07 log gf proved to be consistent with those from other Pr lines. 

We calculated solar photospheric synthetic spectra for all the Pr II lines of Table [6j 
We found, as have the previous studies cited above, that there are few useful solar Pr 
abundance indicators. Our final value was based on five lines (Table [6]). We show the 
synthetic/observed photospheric spectrum matches for four of these lines in left-hand panels 
(a), (c), (e), and (g) of Figure HJ contrasting their appearance in right-hand panels (b), (d), 
(f), and (h) for CS 31082-001. We do not include the 5219.1 A line in Figure 1 because 
it was too weak in the spectrum of CS 31082-001 to analyze in that star. Note that Li07 
transition probabilities were used for the 4222.9, 4510.1, and 5322.8 A lines and IvOl values 
for the 5219.1 and 5259.7 A lines. However, consistent abundances from all five lines were 
derived: the mean value (Table [3]) is log e(Pr) Q = 0.76 ± 0.02 (sigma = 0.04). Our new 
photospheric abundance is in good agreement with the meteoritic and the Biemont et al. 
(1979) photospheric abundances that were quoted above. 

For the r-pro cess-rich low metallicity stars we derived abundances from 10-27 Pr II lines 
(Table [6]). We plot the individual line abundances for these stars and the Sun as functions 
of wavelength in Figure [51 with their summary abundance statistics in the panel legends. In 
each case the line-to-line scatter was small, o ~ 0.06, and we found no significant abundance 
trends with wavelength, excitation energy (the range in this quantity is only ~1 dex), or 

3.2. Dysprosium 

Dy (Z = 66) has seven naturally-occurring isotopes, five of which contribute substan- 
tially to its solar-system abundance: 156 > 158 Dy, <1%, 160 Dy, 2.34%; 161 Dy, 18.91%; 162 Dy, 
25.51%; 163 Dy, 24.9%; and 164 Dy, 28.19% (Lodders 2003). The atomic structure of Dy II is 
complex, leading to a rich spectrum of transitions arising from low-excitation energy levels. 
This species has been well-studied in the laboratory recently, with published transition prob- 
abilities by Kusz (1992), Biemont & Lowe (1993), and Wickliffe, Lawler, and Nave (2000). 
The Wickliffe et al. study contains a detailed comparison of their transition probabilities 
with those of Kusz and Biemont & Lowe (as well as earlier investigations), and will not be 
repeated here. 

We adopted the Wickliffe et al. (2000) log gf values, as in our earlier analyses of the 
r-pro cess- rich stars. Those studies (e.g., Ivans et al. 2006 for HD 221170, and Sneden et al. 
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2003 for CS 22892-052) performed extensive searches for promising Dy II lines. However, the 
Dy abundances reported in those papers were derived from both EW matches and synthetic 
spectrum calculations. Therefore, to be internally consistent in our new analyses, we began 
afresh with new solar Dy II identifications and new synthesis line lists for each chosen feature. 
In principle Dy II lines should have both isotopic wavelength splitting and (for 161 > 163 Dy) 
hyperfine substructure. We inspected the profiles of many of the strongest lines appearing in 
National Solar Observatory (NSO) Fourier Transform Spectrometer (FTS) laboratory Dy II 
spectra. Some line substructure is present in each line. However, the components that are 
shifted away from the line centers are always very weak (<10% of central intensities), and 
the full widths near profile baselines are ~0.05 A. For all lines, FWHM ~ 0.02 A in the lab 
spectra. These widths are substantially smaller than the measured solar and stellar spectrum 
line widths. Therefore we treated all Dy II lines as single features. 

There are many candidate lines, as indicated by their relative strength values shown in 
panel (b) of Figure [TJ Solar Dy abundances could be determined from 13 of these transitions. 
The resulting mean photospheric abundance (Table [3]) is log e(Dy) Q = +1.13 ± 0.02 (a = 
0.06). This value is in excellent agreement with the meteoritic abundance, log e(Dy) met = 
+1.13 ± 0.04 and with the Kusz (1992) photospheric abundance, log e(Dy) = +1.14 ± 0.08. 
It is also in reasonable accord with the Biemont & Lowe (1993) value, log e(Dy) Q = +1.20 ± 0.06. 

Synthetic spectra of 24-35 lines were used in the Dy abundance derivations for the 
r-pro cess- rich low-metallicity giants (Table [7]). The analyses were straightforward, as many 
Dy II lines in each star's spectrum were strong and unblended. This led to very well- 
determined mean abundances (Tables |4] and [5]). 

3.3. Thulium 

Tm [Z = 69) has one naturally-occurring isotope, 169 Tm. This element is one of the 
least abundant of the REs: log e(Tm) mei = 0.11 ± 0.06 (Lodders 2003). Therefore Tm II 
transitions in solar and stellar spectra are weak, and relatively few can be employed in 
abundance analyses. Moore et al. (1966) list only 10 Tm II identifications in their solar line 
compendium; all of these lie at wavelengths A < 4300 A. 

We considered the 146 Tm II lines investigated by Wickliffe & Lawler (1997). That 
study reported laboratory experimental transition probabilities derived from their branching 
fractions and the radiative lifetimes of Anderson, Den Hartog, and Lawler (1996). The 
relative strengths of these lines are displayed in panel (c) of Figure [TJ Inspection of this plot 
suggests that few detectable Tm II lines will be found redward of 4000 A, in accord with the 
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Moore et al. (1966) identifications. 

As in the case of Pr ( §3.11) . we began our search for suitable Tm II transitions with 
CS 31082-001, since they should stand out most clearly among the weaker Fe-peak contam- 
inants in this star's spectrum. Only nine lines were sufficiently strong and unblended to 
warrant further investigation. We computed synthetic spectra for each of these candidate 
features. Although Tm is an odd-Z, odd- A atom with a non-zero nuclear spin (/ = |), 
inspection of the chosen Tm II lines in very high-resolution NSO FTS spectra showed that 
hyperfine splitting is very small, and could be safely ignored in the calculations. 

Our synthetic spectra of Tm II lines for the solar photosphere showed that only three 
of them could be used for abundance analysis. The synthetic/observed spectrum matches 
for these lines in the solar photosphere are displayed in Figure along with those for 
CS 31082-001. It is clear that each of these lines is weak and blended in the photospheric 
spectrum, while being much stronger and cleaner in the r-process-rich low metallicity giant 
star. 

These lines and their photospheric abundances are listed in Table [8J We derive a formal 
mean abundance (Table [3]) of log e(Tm) = +0.14 ± 0.02 (a = 0.04). Caution obviously is 
warranted here. Probably the a value is a truer estimate of the abundance uncertainty than 
the standard deviation of the mean. However, this photospheric abundance is in reasonable 
agreement with the meteoritic value, log e(Tm) mel = +0.11 ± 0.06. 

More Tm II features could be employed in the abundance determinations for the r-process- 
rich low- metallicity giants (Table [8]). Their mean values (Tables H] and [5]) were based on 5-7 
lines per star. For stars analyzed previously by our group, the new Tm abundances agree 
with the published values to within the uncertainty estimates. The Tm abundance for 
CS 31082-OOlwill be discussed along with this star's other REs in §4H 



3.4. Ytterbium 

Yb (Z = 70) has seven naturally-occurring isotopes, six of which are major components 
of its solar-system abundance: 168 Yb, <1%; 170 Yb, 3.04%; m Yb, 14.28%; 172 Yb, 21.83%; 
173 Yb, 16.13%; 174 Yb, 31.83%; and 176 Yb, 12.76% (Lodders 2003). The atomic structure of 
Yb II is similar to that of Ba II, with a 2 S ground state and first excited state more than 
2.5 eV above the ground state. Therefore this species has very strong resonance lines at 
3289.4 and 3694.2 A as the only obvious spectral signatures of this element. All other Yb II 
lines are expected to be extremely weak. 
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The Yb II resonance lines have complex hyperfine and isotopic substructures that 
broaden their absorption profiles by 0.06 A and must be included in synthetic spectrum 
computations. In the Appendix we discuss the literature sources for the resonance lines and 
tabulate their substructures in a form useful for stellar spectroscopists. Moore et al. (1966) 
identified major Fe I, Fe II, and V II contaminants to the 3289.4 A line, and our synthetic 
spectra confirmed that Yb II is a small contributor to the total feature. From our synthetic 
spectra of the 3694.2 A line we derived log e(Yb) = +0.86 ± 0.10 (Table [3]), in reason- 
able agreement with log e(Yb) me t = +0.94 ± 0.03. The large uncertainty attached to our 
photospheric abundance arises from a variety of sources: (a) reliance on a single Yb II line; 
(b) its large absorption strength, which increases the dependence on adopted microturbulent 
velocity; (c) the contaminating presence of the strong Fe I 3694.0 A line; and (b) closeness 
of this spectral region to the Balmer discontinuity. 

We then synthesized the 3289 and 3694 A lines in the stellar spectra. However, these are 
r-pro cess- rich stars, and the isotopic mix in a pure r-process nucleosynthetic mix is different 
than that of the solar system (r-process and s-process) combination. For our computations 
we adopted (see Sneden et al. 2008): ™>™Yb, 0.0%; 171 Yb, 17.8%; 172 Yb, 22.1%; 173 Yb, 
19.0%; 174 Yb, 22.7%; and 176 Yb, 18.4%. 

The Yb contribution to the 3289 A feature is very large in the r-process-rich stars. In 
the most favorable case, CS 31082-001, Yb accounts for roughly 75% of the total blend. 
Unfortunately, the contributions of the contaminants (mostly V II) cannot be assessed accu- 
rately enough for this line to be a reliable Yb abundance indicator. The synthetic/observed 
spectral matches of the 3694 A line provide the new Yb abundances listed in Tables H] and 
[5j These values are consistent with the ones reported in the original papers on these stars. 
However, while the Yb II absorption dominates that of the possible metal-line contaminants, 
the Balmer lines in this spectral region are substantially stronger in these low-pressure giant 
stars than they are in the solar photospheric spectrum. In particular, H I lines at 3691.6 
and 3697.2 A significantly depress the local continuum at the Yb II wavelength. Caution is 
warranted in the interpretation of these Yb abundances. 



3.5. Lutetium 

Lu (Z = 71), has two naturally-occurring isotopes: 175 Lu, 97.416%; and 176 Lu, 2.584% 
(Lodders 2003). It is the least abundant RE: log e(LvL) met = 0.09 ± 0.06 (Lodders). Lu II has 
a relatively simple structure, with a X S ground state. It has no other very low-energy states; 
the first excited level lies 1.5 eV above the ground state. This ion with only two valence 
electrons has relatively few strong lines in the visible and near UV connected to low E.P. 
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levels, although most of the prominent lines have well-determined experimental transition 
probabilities. 

We considered only the Lu II transitions of Quinet et al. (1999), using their experimental 
branching fractions and lifetime measurements by Fedchak et al. (2000) to determine Lu II 
transition probabilities. These are are listed, along with wavelengths and excitation energies, 
in Table 12 of Lawler et al. (2009). The combination of a small solar-system Lu abundance 
and the (unfavorable) atomic parameters produces very small relative strength factors for 
these lines, as shown in panel (d) of Figure [U No line even rises to our defined "weak-line" 
threshold of usefulness. Moore et al. (1966) lists only 3077.6, 3397.1, and 3472.5 A Lu II 
identifications in their solar line compendium, and all of these lines appear to be blended. 

We made a fresh search for detectable lines of Lu II, and succeeded mainly in confirming 
the results of a previous investigation by Bord, Cowley, & Mirijanian (1998). Those authors 
argued that all of the lines identified by Moore et al. (1966) are unsuitable for solar Lu 
abundance work. They quickly dismissed the 3077.6 and 3472.5 A lines and performed an 
extended analysis of 3397.1 A. Synthetic spectrum computations around this feature (see 
their Figures 2 and 3) convinced them that molecular NH dominates the absorption at the 
Lu II wavelength. Our own trials produced the same outcome. 

Bord et al. (1998) detected Lu II 6221.9 A in the Delbouille et al. (1973) photospheric 
spectrum. This line is extremely weak, EW ~ 1 mA, and its hyperfine substructure spreads 
the absorption over ~0.5 A. The complex absorption profile of this line (see their Figure 4) 
actually increases one's confidence in its identification in the photospheric spectrum. Bord 
et al. reported log e(Lu) = +0.06 with an estimated ±0.10 uncertainty from this line. 

We repeated their analysis, using the hyperfine substructure pattern given in Table 13 
of Lawler et al. (2009), and derived log e(Lu) Q = +0.12 ± 0.08 (Table [3]), where the error 
reflects uncertainties in matching synthetic and observed feature profiles. This photospheric 
abundance is consistent with the Bord et al. (1998) value and with the meteoritic abundance 
quoted above, given the uncertainties attached to each of these estimates. Our lack of success 
in identifying other Lu abundance indicators in the solar photospheric spectrum suggests that 
prospects are poor for reducing its error bar substantially in the future. 

We also attempted to study the 3397 and 6621 A lines in our sample of r-process-rich 
low metallicity giants. Absorption by Lu II at 3397.1 A is certainly present in the spectra 
of at least CS 22892-052 and CS 31082-001. Unfortunately, the lower resolutions of our 
stellar spectra compared to that of the solar spectrum creates more severe blending of the 
Lu transition with neighboring lines, and NH contamination of the total feature still creates 
substantial abundance ambiguities. The 6221.9 A line should be present, albeit very weak, 



13 



in these stars. However, our spectra (when they extend to this wavelength range) lack the 
S/N to allow meaningful detections. We therefore cannot report Lu abundances for these 
r-pro cess- rich stars. 

4. RARE EARTH ABUNDANCE DISTRIBUTIONS IN THE SUN AND 

R-PROCESS-RICH STARS 

4.1. The Sun and Solar System 

With new analyses of Pr, Dy, Tm, Yb, and Lu we now have determined abundances 
for the entire suite of REs in the solar photosphere. In Table |3] we merge the results of this 
and our previous papers. Missing from the list is of course Pm (Z = 61), whose longest- 
lived isotope, 145 Pm, is only 17.7 years (Magill, Pfennig, and Galy 2006). We also chose 
not to include a photospheric value for Ba, whose few transitions are so strong that their 
solar absorptions cannot be reliably modeled in the sort of standard photospheric abundance 
analysis that we have performed. 

The photospheric abundance uncertainties quoted in Table [3] are combinations of in- 
ternal "scatter" factors (mainly continuum placement, observed/synthetic matching, and 
line blending problems) and external "scale" factors (predominantly solar model atmosphere 
choices). These issues are discussed Lawler et al. (2009) and in previous papers of this series. 
We remind the reader that our abundance computations have been performed with the tra- 
ditional assumptions of LTE and ID static atmosphere geometry. Very little has been done 
to date to explore the effects of these computational limitations for RE species in the solar 
atmosphere. Mashonkina & Gehren (2000) have performed non-LTE abundance analyses of 
Ba and Eu, but their photospheric abundances are not substantially different from LTE re- 
sults. There have been efforts to model the solar spectrum with more realistic 3-dimensional 
(3D) hydrodynamic models; see the summary in Grevesse, Asplund, & Sauval (2007), and 
references therein. These studies so far have reported new solar abundances only for the 
lighter elements (CNO, Na— Ca, and Fe). Generally the 3D non-LTE line modeling efforts 
yield lower abundances: comparing the photospheric values in Grevesse et al. to those of 
the older standard compilation of Anders & Grevesse (1989), < <51og e > = —0.12 ± 0.03 
(a = 0.09, for 11 elements that can be studied with photospheric spectra). We thus expect 
that any RE abundance shifts with 3D modeling would be similar from element to element, 
leaving their abundance ratios essentially unchanged. Future studies to explore these effects 
in detail will be welcome. 

In Figure [7] we compare RE photospheric abundances to their meteoritic values. In 
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the top panel the "OLD" values are best estimates by Anders & Grevesse (1989). While 
the average agreement is good, significant discrepancies between individual abundances are 
evident, particularly at the low-abundance end. Formally, a simple mean is < loge _^G89 ~ 
loge met _ AG89 > = 0.00 ± 0.06 (a = 0.22). In the bottom panel, the "NEW" meteoritic 
abundances (Lodders 2003) are correlated with our "NEW" photospheric ones (Table[3]). The 
data sources are denoted by different symbols in the figure: red open circles for photospheric 
abundances newly determined here and in Lawler et al. (2009) for which Wisconsin-group 
lab data have been used; black filled circles for abundances reported in our previous papers; 
and blue open triangles for two elements with transition probability data adopted from 
other literature sources. Clearly the agreement is excellent: for 15 elements the formal mean 
difference is < loge Q — loge met > = 0.01 ± 0.01 (er = 0.05). No trends are discernible with 
the source of atomic data, or the abundance levels (as shown in the figure), or the number of 
lines that contribute to the photospheric abundances (Tabled]). With the possible exception 
of Hf (discussed in Lawler et al. 2007 and in §5]), and with repeated cautions about the 
photospheric abundances deduced from only one or two transitions, the two primary sources 
of primordial Solar-System abundances appear to be in complete accord. 



4.2. The r-process-Rich Low Metallicity Giant Stars 

Rare-earth abundances for the five r-process-rich stars from this and our previous papers 
are collected in Tables H] and [5j For all stars the Pr, Dy, Tm, and Yb abundances are, of 
course, newly determined in this paper. We chose also to redo all the Ba abundances via 
new synthetic spectrum calculations, to ensure that these were determined in a consistent 
manner. We also performed new analyses for selected elements in individual stars (e.g., Tb 
in HD 115444) when the original papers either did not report abundance values or did so 
with now-outdated atomic data. 

Of particular interest is the very n-capture-enhanced star CS 31082-001, which is a 
recent addition to our r-process-rich star list. This star gained notoriety as the first r-process- 
rich star with a convincing detection of U, a long-lived radioactive element of great interest 
to cosmochronology (Cayrel et al. 2001). The first and most complete study of this star was 
published by Hill et al. (2002). The mean difference between our RE abundances for this star 
and theirs is < logenui — loge us > = —0.05 ± 0.03 (cr = 0.10, 12 elements in common). We 
also compared our CS 31082-001 abundances with those of Honda et al. (2004), with similar 
results: < log enonda — loge us > = +0.07 ± 0.03 (a = 0.09, 12 elements in common). The 
mean offsets are very small, and reflect minor differences in model atmospheres, observed 
spectra, analytical techniques, and atomic data choices. The element-to-element scatters are 
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also reasonable, given the use of many more transitions in our study (a total of 342, Table H]) 
compared to 95 in Hill et al. and 49 in Honda et al. Note that some portion of the a's in 
these comparisons arises because the Tb abundance differences are offset by ~0.2 dex from 
the mean differences (we derive larger values). Investigation of this one anomaly is beyond 
the scope of this work. 

The abundance standard deviations of samples (a) and of means that are given in 
Tables H] and refer to internal (measurement scatter) errors only. To investigate scale 
uncertainties, we determined the abundance sensitivities of eight RE elements to changes 
in model parameters (T e ff, log g, [M/H], £ t )> to changes in the adopted model atmosphere 
grid, and to changes to line computations to better account for continuum scattering opac- 
ities. In Table M we summarize the results of these exercises. We began with a "baseline" 
model atmosphere from the Kurucz (1998) grid with parameters T e g = 4750 K, log g = 1.5, 
[M/H] = -2.5, and £ t — 2.0. Such a model is similar to the ones adopted for the r-process-rich 
giants (Table [2]). We derived abundances with this model for 1-4 typical transitions each of 
the elements for the program star CS 31082-001. Full account was taken of hyperfme and 
isotopic substructure for La, Pr, Eu, and Yb. We then repeated the abundance derivations 
for models with parameters varied as indicated in Table [HJ including a trial using a model 
with baseline parameters taken from the new MARCS (Gustafsson et al. 2008) gri d@ The 
inclusion of scattering in computations of continuum source functions, a new feature in our 
analysis code, is described in Sobeck et al. (2008) 

The Table M quantities are differences between abundances of the individual models 
and those of the baseline model. The uncertainties in stellar model parameters given in the 
original r-process-star papers are typically ±150 K in T c fj, ±0.3 in log g, ±0.2 km s _1 in £ t , 
and ±0.2 in [M/H] metallicity. Application of these uncertainties to the model parameter 
dependences of Table M suggests that [M/H] and £ t choices are not important abundance 
error factors. Temperature and gravity values obviously play larger roles. However, while 
the absolute abundances of individual elements change with different T e g and log g choices, 
the relative abundances generally do not; in most cases, all RE abundances move in lock step. 
Assuming here that the atmosphere parameter uncertainties are uncorrelated, we estimate 
total abundance uncertainties for each RE element to be ~0. 15— 0.20, but the abundance 
ratios have uncertainties of ~0. 01— 0.05 (the exception is Yb, represented by only one very 
strong line in the UV spectral region; see §3.4p . More detailed computations that consider 
departures from LTE among RE first ions in the atmospheres of very metal-poor giant stars 
should be undertaken in the future. Some first steps in this direction have been undertaken 
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for Ba and Eu by Mashonkina et al. (2008), but such calculations will need to be repeated 
for many REs to understand the magnitude of corrections to the abundances reported here. 

5. DISCUSSION 

We illustrate the RE abundances for BD+17 3248, CS 22892-052, CS 31082-001, HD 115444 
and HD 221170 in Figures [8] and [9j For each star the abundances have been normalized at 
Eu, a predominantly r-process element. In Figure [8] these relative abundances are shown in 
comparison to the Solar System r-process-only predictions from Arlandini et al. (1999) and 
Simmerer et al. (2004). We note first the excellent star-to-star (relative abundance) agree- 
ment. Early RE abundance distributions of n-capture-rich metal-poor stars indicated large 
star-to-star scatter for a number of individual elements (e.g., Luck & Bond 1985, Gilroy et al. 
1988). The combination of substantially better S/N and resolution of the stellar spectra and 
the experimental initiatives of this series of papers has dramatically reduced that scatter - 
all the RE elements are now in very good (relative) agreement for these five halo stars. 

Figure [8] also uses solid lines to illustrate the solar-system r-process-only meteoritic 
abundances determined by Simmerer et al. (2004) and Arlandini et al. (1999). In both 
cases, these values were computed by subtracting the s-pro cess-only abundances from the 
total elemental abundances. The "classical" method (Simmerer et al.) matches smooth aN s 
curves to those isotopes of n-capture elements whose production is essentially all due to the 
s-process, and infers from those empirical curves the s-process amounts of elements that can 
be produced by both the r-process and s-process. The solar system r-process abundances are 
then just the residuals between total elemental and s-process amounts. The "stellar" method 
(Arlandini et al.) uses theoretical models of s-process nucleosynthesis instead of empirical 
s-process abundance curves, and again infers the r-process amounts by subtraction. 

Our stellar abundances compare very well with the relative solar system r-process dis- 
tributions. In the past we and other investigators have found overall agreement, but on a 
more approximate scale. The new abundance determinations shown in Figure [H] tighten the 
comparison, with deviations between the stellar and solar system r-process curves of typi- 
cally less than 0.1 dex - probably the practical limit of what is currently possible. These 
abundance comparisons strongly support many other studies (see Sneden et al. 2008, and 
references therein) arguing that essentially the same process was responsible for the forma- 
tion of all of the r-process contributions to these elements early in the history of the Galaxy 
in the element progenitor stars to the presently-observed r-process-rich halo stars. 

Despite this general level of elemental abundance consistency, there are some interesting 
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deviations. In particular, the two solar system r-process predictions differ by about 0.1 dex 
for the elements Ce and Nd (Table [3]). In both cases the stellar model predictions from 
Arlandini et al. (1999) give a better fit to the stellar abundance data than do the stan- 
dard model predictions from Simmerer et al. (2004). This suggests that the Arlandini et al. 
r-process distribution might be superior for such abundance comparisons. This has been 
noted previously by others (e.g., Roederer et al. 2008) for isotopic studies. 

There is also still some star-to-star scatter particularly at Ba, with several stellar ele- 
mental abundances appearing somewhat higher than the solar system r-process curves. This 
can be seen more clearly in Figure [9l where we illustrate the difference between the relative 
(scaled to Eu) stellar RE and the scaled solar system r-process abundances (Arlandini et al. 
1999) in the five r-pro cess-rich stars. While most of the individual elemental abundance 
data lie close to the dotted line (indicating perfect agreement with the solar r-process), Ba 
and Yb have significant star-to-star scatter. But both elements have inherent observational 
problems, as they are represented by only a few very strong transitions that have multi- 
ple isotopic components whose relative abundances are sensitive to the relative r-/s-process 
dominance (recall the Yb discussion in I3.4D . Abundance determinations for Yb and Ba are 
less reliable than those of most other RE elements, and should be treated with caution. 

We also note that for BD+17 3248 the RE abundances relative to Eu appear to be 
somewhat higher their values in the other stars, particularly for the predominantly s-process 
elements Ba and La. BD+17 3248 has a metallicity of [Fe/H] ~ -2.1 (Cowan et al. 2002), 
so this star might be showing the signs of the onset of Galactic s-processing, which occurs 
at approximately that metallicity (Burris et al. 2000). On the other hand HD 221170 with 
a similar metallicity (Ivans et al. 2006) does not seem to show the same deviations for the 
s-process elements, and thus the deviations for BD+17 3248 may be specific to that star. 

We examine whether there is any correlation between the deviation of the stellar abun- 
dances from the solar system r-process values and the s-process percentage of those elements 
in solar system material (from Simmerer et al. 2004) in Figure [TOj It is clear that there is lit- 
tle if any secular trend with the abundance differences with increasing solar-system s-process 
abundance percentage. This lack of correlation was also found specifically for the element 
Ce by Lawler et al. (2009). 

To get a clearer sense of the overall abundance agreement with the solar- system r-process 
abundances, we show in Figure [TT1 the arithmetic averages of the elemental abundance offsets 
(from Figure [TO]) for the five stars, again as a function of s-process percentage. Obviously 
these average offsets with respect to the solar-system r-process values are very small. Includ- 
ing all elements the mean of the average offsets is loge = 0.05 (a = 0.05). Previously Lawler 
et al. (2007) had found that the observed average stellar abundance ratio of Hf/Eu in a group 
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of metal-poor halo stars is larger than previous estimates of the solar-system r-process-only 
value, suggesting a somewhat larger contribution from the r-process to the production of 
Hf. Our new analysis supports that finding, as the average Hf offset is larger than all of the 
other elemental abundances. If the solar system r-process contribution was larger it would 
drive down the average offset illustrated in Figure [TTJ Ignoring the Hf results, the mean of 
the average offsets for all of the other RE elements is 0.04 (a = 0.03). This is essentially 
a perfect agreement within the limits of our observational and experimental uncertainties, 
as well as the uncertainties (observational and theoretical) associated with the solar system 
r-process-only abundance values. 



6. CONCLUSIONS 

We have determined new abundances of Pr, Dy, Tm, Yb, and Lu for the solar pho- 
tosphere and for five very metal-poor, r-process-rich giant stars. Combining these results 
with those of previous papers in this series (cited in §1), we have now derived very accurate 
solar/stellar abundances for the entire suite of stable RE elements. 

With the single exception of Hf, the solar photospheric abundances agree with solar- 
system meteoritic values perfectly to within the uncertainty estimates of each. Our photo- 
spheric and stellar analyses have emphasized studying as many transitions of each species 
as possible (up to 46 Nd II lines in the Sun, up to 72 Sm II lines in BD+17 3248). The 
line-to-line abundance scatters are always small when the number of available transitions 
is large (typically o < 0.07). This clearly demonstrates the reliability of the RE transition 
probabilities published in this series of papers. We argue that, with proper care in stellar 
analyses, trustworthy abundances of RE elements can be now be determined from spectra 
in which far fewer transitions are available. 

Utilizing the new experimental atomic data we have determined far more precise stellar 
RE elemental abundances in five r-process rich stars. These newly derived values show a 
dramatic decrease in star-to-star elemental abundance scatter - all the RE elements are 
now in very good (relative) agreement for these five halo stars. Furthermore, our newly 
derived values indicate an almost perfect agreement between the average stellar abundances 
and the solar system r-process-only abundances for a wide range of elements in these five 
r-process-rich stars. There is no evidence for significant s-process contamination. The one 
exception appears to be a somewhat higher value of stellar Hf with respect to the solar 
system r-process-only value for this element. This may indicate that further analysis of the 
solar r- and s-process deconvolution for this element might be useful. These results for the 
five r-process-rich halo stars confirm, and strongly support, early studies that indicate that 



-19- 



the r-process was dominant for the n-capture elements early in the history of the Galaxy. 
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tronomico di Padova; the Director and staff are thanked for their hospitality and financial 
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solar abundance analyses of the present and previous papers of this series have greatly ben- 
efited from the availability of the photospheric spectrum in the BASS2000 Solar Survey 
Archive maintained by l'Observatoire de Paris. This work has been supported by the Na- 
tional Science Foundation through grants AST 05-06324 to JEL and EDH, AST 06-07708 
to CS, and AST 07-07447 to JJC. 



-20- 



APPENDIX 

There have been numerous experimental studies of hyperfine structure (hfs) in Pr II. 
We have reviewed the literature for measurements on the upper and lower levels of lines 
useful, or potentially useful, for elemental abundance studies. Six publications are relevant, 
as indicated in Table [101 One sees generally good agreement among measured values of the 
hfs A constants. Only a few, not very accurate, measurements of the hfs B constants have 
been reported. Since the electric quadrupole interaction (B constants) has a much smaller 
effect on the line component pattern than the magnetic dipole interaction (A constant), it 
is often neglected and will be neglected here. 

One of the best and fairly extensive set of measurements of Pr II hfs A constants is that 
by Rivest et al. (2002) using laser induced fluorescence. We adopted their measurements, if 
available, to compute the complete hfs line component patterns that are given in Table [TT1 
For levels which were not studied by Rivest et al., we used hfs A constants from Ginibre 
(1989). IvOl improved some Pr II energy levels using FTS data. The center-of-gravity 
wavenumbers in Table [TT] are from the IvOl energy levels in every case where an improved 
energy was reported for both the upper and lower level of the line. For other lines the 
center-of gravity wavenumbers are from the NIST energy levels (Martin et al. 1978), because 
it is probably not a good idea to mix energy levels from two sources. Center-of-gravity air 
wavelengths were computed from wavenumbers using the standard index of air (Edlen 1953). 

For Yb II we used the isotopic and hyperfine data of Martensson-Pendrill, Gough, & 
Hannaford (1994). We adopted the transition probabilities of Biemont et al. (1998) renor- 
malized to the lifetime results of Pinnington, Rieger, & Kernahan (1997): log gfmm = +0.02, 
and log gf = —0.30. These values are close to those derived from Biemont et al. (2002), 
as given in the D.R.E.A.M. databas^: log (7/3289 = —0.05, and log (7/3694 = —0.32. Combin- 
ing the transition probabilities, hyperfine and isotopic substructures, and the solar isotopic 
breakdown given in 13.41 yields complete transition structures for these two Yb II lines; these 
are listed in Table [121 



http: / / w3.umh.ac.be/~astro/dream.shtml 



- 21 - 



REFERENCES 

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 

Anderson, H. M., Den Hartog, E. A., & Lawler, J. E. 1996, J. Opt. Soc. America B, 13, 2382 

Arlandini, C, Kappeler, R, Wisshak, K., Gallino, R., Lugaro, M., Busso, M., & Straniero, 
O. 1999, ApJ, 525, 886 

Bernstein, R., Shectman, S. A., Gunnels, S. M., Mochnacki, S., & Athey, A. E. 2003, 
Proc. SPIE, 4841, 1694 

Biemont, E., Dutrieux, J.-R, Martin, I., & Quinet, P. 1998, Journal of Physics B Atomic 
Molecular Physics, 31, 3321 

Biemont, E., Grevesse, N., & Hauge, O. 1979, Sol. Phys., 61, 17 

Biemont, E., Lefebvre, P.-H., Quinet, P., Svanberg, S., & Xu, H. L. 2003, European Physical 
Journal D, 27, 33 

Biemont, E., & Lowe, R. M. 1993, A&A, 273, 665 

Biemont, E., Quinet, P., Dai, Z., Zhankui, J., Zhiguo, Z., Xu, H., & Svanberg, S. 2002, 
Journal of Physics B Atomic Molecular Physics, 35, 4743 

Bord, D. J., Cowley, C. R., & Mirijanian, D. 1998, Sol. Phys., 178, 221 

Burris, D. L., Pilachowski, C. A., Armandroff, T. E., Sneden, C, Cowan, J. J., & Roe, H. 
2000, ApJ, 544, 302 

Cayrel, R., et al. 2001, Nature, 409, 691 

Cowan, J. J., et al. 2005, ApJ, 627, 238 

Cowan, J. J., et al. 2002, ApJ, 572, 861 

Delbouille, L., Roland, G., & Neven, L. 1973, Atlas photometrique du spectre solaire de 
A3000 a A10000, Liege: Universite de Liege, Institut d'Astrophysique 

Den Hartog, E. A., Herd, M. T., Lawler, J. E., Sneden, C, Cowan, J. J., & Beers, T. C. 
2005, ApJ, 619, 639 

Den Hartog, E. A., Lawler, J. E., Sneden, C, & Cowan, J. J. 2003, ApJS, 148, 543 
Den Hartog, E. A., Lawler, J. E., Sneden, C, & Cowan, J. J. 2006, ApJS, 167, 292 



-22 - 



Edlen, B. 1953, Journal of the Optical Society of America, 43, 339 

Fedchak, J. A., Den Hartog, E. A., Lawler, J. E., Palmeri, P., Quinet, P., & Biemont, E. 
2000, ApJ, 542, 1109 

Fitzpatrick, M. J., & Sneden, C. 1987, BAAS, 19, 1129 

Gilroy, K. K., Sneden, C, Pilachowski, C. A., & Cowan, J. J. 1988, ApJ, 327, 298 
Ginibre, A. 1989, Phys. Scr, 39, 694 

Grevesse, N., Asplund, M., & Sauval, A. J. 2007, Space Science Reviews, 130, 105 

Gustafsson, B., Edvardsson, B., Eriksson, K., J0rgensen, U. G., Nordlund, A., & Plez, B. 
2008, A&A, 486, 951 

Heifer, H. L., Wallerstein, G., & Greenstein, J. L. 1959, ApJ, 129, 700 

Hill, V., et al. 2002, A&A, 387, 560 

Holweger, H., & Miiller, E. A. 1974, Sol. Phys., 39, 19 

Honda, S., Aoki, W., Kajino, T., Ando, H., Beers, T. C, Izumiura, H., Sadakane, K., & 
Takada-Hidai, M. 2004, ApJ, 607, 474 

Ivans, I. L, Simmerer, J., Sneden, C, Lawler, J. E., Cowan, J. J., Gallino, R., & Bisterzo, 
S. 2006, ApJ, 645, 613 

Ivarsson, S., Litzen, U., & Wahlgren, G. M. 2001, Phys. Scr, 64, 455 (IvOl) 

Ivarsson, S., Wahlgren, G. M., & Ludwig, H.-G. 2003, BAAS, 35, 1421 

Kelson, D. D. 2003, PASP, 115, 688 

Kelson, D. D., Illingworth, G. D., van Dokkum, P. C, & Franx, M. 2000, ApJ, 531, 159 

Kurucz, R. L. 1998, Fundamental Stellar Properties, IAU Symposium 189, (Dordrecht: 
Kluwer), 217 

Kusz, J. 1992, A&AS, 92, 517 

Lawler, J. E., Bonvallet, G., & Sneden, C. 2001a, ApJ, 556, 452 

Lawler, J. E., Den Hartog, E. A., Sneden, C, & Cowan, J. J. 2006, ApJS, 162, 227 



-23- 



Lawler, J. E., Den Hartog, E. A., Sneden, C, & Cowan, J. J. 2008a, Canadian J. Phys., 86, 
1033 

Lawler, J. E., Hartog, E. A. D., Labby, Z. E., Sneden, C, Cowan, J. J., & Ivans, I. I. 2007, 
ApJS, 169, 12 

Lawler, J. E., Sneden, C, k Cowan, J. J. 2004, ApJ, 604, 850 

Lawler, J. E., Sneden, C, Cowan, J. J., Wyart, J.-F., Ivans, I. I., Sobeck, J. S., Stockett, 
M. H., & Den Hartog, E. A. 2008b, ApJ, 178, 71 

Lawler, J. E., Wickliffe, M. E., Cowley, C. R., & Sneden, C. 2001b, ApJS, 137, 341 

Lawler, J. E., Wickliffe, M. E., Den Hartog, E. A., & Sneden, C. 2001c, ApJ, 563, 1075 

Lawler, J. E., Sneden, C, Cowan, J. J., Ivans, I. I., & Den Hartog, E. A. 2009, ApJS, in 
press 

Li, M., Ma, H., Chen, M., Chen, Z., Lu, F., Tang, J., & Yang, F. 2000a, Hyperfine Interac- 
tions, 128, 417 

Li, M., Ma, H., Chen, M., Lu, F., Tang, J., & Yang, F. 2000b, Phys. Rev. A, 62, 052504 

Li, R., Chatelain, R., Holt, R. A., Rehse, S. J., Rosner, S. D., & Scholl, T. J. 2007, Phys. Scr, 
76, 577 (Li07) 

Lodders, K. 2003, ApJ, 591, 1220 

Luck, R. E., & Bond, H. E. 1985, ApJ, 292, 559 

Ma, H., Chen, M., Chen, Z., Shi, W., Lu, F., Tang, J., & Yang, F. 1999, Journal of Physics 
B Atomic Molecular Physics, 32, 1345 

Magill, J., Pfennig, C, & Galy, J. 2006, Karlsruher Nuklidkarte, Chart of the Nuclides. 

Malcheva, C, Blagoev, K., Mayo, R., Ortiz, M., Xu, H. L., Svanberg, S., Quinet, P., & 
Biemont, E. 2006, MNRAS, 367, 754 

Martensson-Pendrill, A.-M., Gough, D. S., & Hannaford, P. 1994, Phys. Rev. A, 49, 3351 

Martin, W. C, Zalubas, R., & Hagan, L. 1978, Atomic Energy Levels - The Rare-Earth El- 
ements, NSRDS-NBS, Washington: National Bureau of Standards, U.S. Department 
of Commerce, 1978 

Mashonkina, L., & Gehren, T. 2000, A&A, 364, 249 



-24- 



Mashonkina, L., et al. 2008, A&A, 478, 52 

Moore, C. E., Minnaert, M. G. J., & Houtgast, J. 1966, The solar spectrum 2935 A to 
8770 A, Nat. Bur. Stds. Monograph, Washington: US Government Printing Office 

O'Brien, S., Dababneh, S., Heil, M., Kappeler, F., Plag, R., Reifarth, R., Gallino, R., & 
Pignatari, M. 2003, Phys. Rev. C, 68, 035801 

Pinnington, E. H., Rieger, G., & Kernahan, J. A. 1997, Phys. Rev. A, 56, 2421 

Quinet, P., Palmeri, P., Biemont, E., McCurdy, M. M., Rieger, G., Pinnington, E. H., 
Wickliffe, M. E., & Lawler, J. E. 1999, MNRAS, 307, 934 

Rivest, R. C, , Izawa, M. R., Rosner, S. D., Scholl, T. J., Wu, G., & Holt, R. A., 2002, 
Canadian Journal of Physics, 80, 557 

Roederer, I. U., Lawler, J. E., Sneden, C, Cowan, J. J., Sobeck, J. S., & Pilachowski, C. A. 
2008, ApJ, 675, 723 

Simmerer, J., Sneden, C, Cowan, J. J., Collier, J., Woolf, V. M., & Lawler, J. E. 2004, ApJ, 
617, 1091 

Sneden, C. 1973, ApJ, 184, 839 

Sneden, C, et al. 2003, ApJ, 591, 936 

Sneden, C, Cowan, J. J., & Gallino, R. 2008, ARA&A, 46, 241 

Sobeck, J. S., Kraft, R. P., and Sneden, C. 2008, A J, to be submitted 

Sobeck, J. S., Lawler, J. E., & Sneden, C. 2007, ApJ, 667, 1267 

Tull, R. C, MacQueen, P. J., Sneden, C, & Lambert, D. L. 1995, PASP, 107, 251 

Vogt, S. S., et al. 1994, Proc. SPIE, 2198, 362 

Westin, J., Sneden, C, Gustafsson, B., k Cowan, J. J. 2000, ApJ, 530, 783 

Wickliffe, M. E., & Lawler, J. E. 1997, J. Opt. Soc. America B, 14, 737 

Wickliffe, M. E., Lawler, J. E., & Nave, G. 2000, Journal of Quantitative Spectroscopy and 
Radiative Transfer, 66, 363 



This preprint was prepared with the AAS IATgX macros v5.2. 



-25 - 



2 

1 



■1 

-2 

-3 
2 

1 



-1 

■2 

■3 



Pr 


II 




strong lines 


• 


i | detection _ _ ■ 
' limit 


« •• • > • • . • 

• 

• * 


. • / • . • 

• * 

% « * ■. • ♦ 

f : • . 
• • • • • 

• • 


" (a) 

i 


• • • 

i 


4000 


6000 



Tm II 



Dy II 



' . V • • 4 



• % •• • 



- •-• *<$. : V. 

• _ I ■ • 
• • • • • 

• • :*. .• 
- % 



(b) 



_L 



_L 



4000 



6000 



Lu II 



» ♦ 



(c) 



4000 



(d) 



6000 4000 
Wavelength (A) 



6000 



Fig. 1. — Relative strength plots for Pr II, Dy II, Tm II, and Lu II lines. Dashed lines in each box denote 
the approximate lower strength limit for strong lines, and dotted lines denote the lower limit for detectable 
lines, as defined in the text. For these plots the wavelength range has been restricted to A > 3000 A (the 
cutoff for ground-based spectra) and A > 7500 A (for lack of detectable lines of these species) . 
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Fig. 2.— Differences between Ivarsson et al. (2001; IvOl) and Li et al. (2007; Li07) Pr II 
log gf values plotted as a function of wavelength. As indicated in the figure legend, red dots 
denote transitions employed in our solar/stellar analyses, and blue x symbols denote other 
lines in common between IvOl and Li07. 




Wavelength (A) 

Fig. 3. — Observed and synthetic spectra in CS 31082-001 of two strong Pr II lines with 
wide hyperfine structure. In each panel, the points represent the observed spectrum. The 
magenta line is the spectrum computed with no contribution from Pr II; the black line is the 
best-fitting synthesis (with the Pr abundance given in Table El and the red and blue lines 
are the syntheses computed with Pr abundances altered by ±0.3 dex from the best value. 
Vertical lines have been drawn at the bottom of each panel to indicate the wavelengths 
and relative strengths (arbitrary overall normalization) of the hyperfine components that 
comprise the Pr II transitions. 
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Fig. 4. — Observed and synthetic spectra of the Sun (left-hand panels (a), (c), (e), (g)) and 
CS 31082-001 (right-hand panels (b), (d), (f), (h)) for the four Pr II lines that contribute to 
the solar abundance estimate. In each panel, the points represent the observed spectrum. 
The magenta line is the spectrum computed with no contribution from Pr II; the black line 
is the best-fitting synthesis (with the Pr abundance given in Table E]); and the red and blue 
lines are the syntheses computed with Pr abundances altered by ±0.3 dex from the best 
value. The solar spectrum is that of Delbouille et al. (1973), but sampled at a wavelength 
step size of 0.01 A for display purposes. 
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Fig. 6. — Observed and synthetic spectra of the Sun (left-hand panels (a), (c), (e)) and CS 31082-001 
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from the best value. The solar spectrum is that of Delbouille et al. (1973), but sampled at a wavelength step 
size of 0.01 A for display purposes. 
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Fig. 7. — Comparisons of solar-system meteoritic and solar photospheric abundances of the 
RE elements. In the top panel, the "OLD" abundances are the recommended values from 
two decades ago (Anders & Grevesse 1989). In the bottom panel, the "NEW" meteoritic 
values are from (Lodders 2003), and the solar photospheric abundances are from in this 
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elements whose abundances are based on other sources for the transition probabilities, and 
black dots for abundances determined in earlier papers of this series. 
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Fig. 8. — Comparison of the newly derived RE abundances in five r-rich halo stars with 
predictions for solar system r-process only abundances from Arlandini et al. (1999) and 
Simmerer et al. (2004). For each star the abundances have been normalized at Eu. 



-33- 



0.50 



0.30 



0.10 



-0.10 



-0.30 



-0.50 




La 

Ce Pr Nd 




Hf 



Yb 



Sm Gd 



► BD+17 3248 
■ CS 22892-052 

• CS 31082-001 

HD1 15444 
AHD221170 

- SS r-process 



55 



65 



75 



Atomic Number 



Fig. 9. — Comparison of the newly derived RE abundances in five r-rich halo stars to the 
solar system r-process only value from Arlandini et al. (1999). For each star the abundances 
have been normalized at Eu. The dotted line indicates a perfect agreement between the 
stellar and solar system r-only values. The error bars are the sigma values listed for each 
star in Tables H] and [H 
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Fig. 10. — Comparison of the newly derived RE abundances in five r-rich halo stars to the 
solar system r-process only value from Arlandini et al. (1999) as a function of percentage of 
the solar system elemental s-process. The dashed line indicates a perfect agreement between 
the stellar and solar system r-only values. For clarity in this figure, a different color has been 
used for each element. 
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Fig. 11. — Averages of the stellar elemental abundance offsets of the five stars with respect to 
the solar system r-process only value from Arlandini et al. (1999) as a function of percentage 
of the solar system elemental s-process. For each star included in the average the abundance 
offsets have been normalized at Eu. The dotted line indicates a perfect agreement between 
the average stellar and solar system r-only values. 
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Table 1. Spectroscopic Observations 



Snprrropra dIi 


A Ran pp 

A 




S/N 


F h 


A c 

"app 

A 


Stars 


Keck I HIRES' 1 


3050-5950 


40,000 


100 


1142 


3500 


CS 31082-001, HD 221170 








150 


1500 


4000 


CS 31082-001, HD 221170 








200 


1778 


4500 


CS 31082-001, HD 221170 




3100-4250 


45,000 


100 


1286 


3500 


CS 22892-052 








150 


1688 


4000 


CS 22892-052 








200 


2000 


4500 


CS 22892-052 




3100-4650 


45,000 


100 


1286 


3500 


BD+17 3248, HD 115444 








150 


1688 


4000 


BD+17 3248, HD 115444 








200 


2000 


4500 


BD+17 3248, HD 115444 


McDonald "2d-coude" c 


3800-9000 


60,000 


100 


1500 


4000 


BD+17 3248, HD 115444 








250 


2500 


6000 


BD+17 3248, HD 115444 




3800-7800 


60,000 


80 


1200 


4000 


HD 221170 








275 


2750 


6000 


HD 221170 




4200-7000 


60,000 


50 


750 


4000 


CS 22892-052 








150 


1500 


6000 


CS 22892-052 


Magellan Clay MIKE f 


3800-4950 


50,000 


100 


1250 


4000 


CS 22892-052 




5050-8000 


38,000 


150 


950 


6000 


CS 22892-052 



a R = X/SX 

h F = (R/X app x (S/N)) 

c X app = the approximate wavelength for calculation of F 

d Vogt ct al. (1994); detailed description at http://www.ucolick.org/~hires/ 

c Tull et al. (1995); detailed description at http://www.as.utexas.edU/mcdonald/facilities/2.7m/cs2.html 
f Bernstein et al. (2003); detailed description at http://www.ucolick.org/~rab/MIKE/usersguide.html 
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Table 2 




ellar Moc 


iel Parameters 


Star 


T eff 


log 9 


[Fc/H] 


v t 


Reference 




K 






km s _1 




BD+17 3248 


5200 


1.80 


-2.10 


1.90 


Cowan ct al. (2002) 


CS 22892-052 


4800 


1.50 


-3.12 


1.95 


Sneden et al. (2003) 


CS 31082-001 


4825 


1.50 


-2.91 


1.90 


Hill et al. (2002) 


HD 115444 


4800 


1.50 


-2.90 


2.00 


Westin et al. (2000) 


HD 221170 


4510 


1.00 


-2.19 


1.80 


Ivans et al. (2006) 
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Table 3. Solar Rare Earth Abundances 



Element 


Z 


log 


5 € me t 




lo 




© 




a 




Rcf° 


log TV/ 


log TV/ 






meteoritic 


















empirical 


stellar 


Ba 


56 


2.19 


± 


0.03 
















1 


-0.0936 


-0.0696 


La 


57 


1.18 


± 


0.06 


1, 


.14 


± 


0. 


.01 


0.03 


14 


2 


-0.9547 


-0.9210 


Cp 


oo 


1.61 


± 


0.02 


i 


. yj l 




n 

u 


m 


u.uu 


4: 


Q 
o 




U .O i oo 


Pr 


59 


0.78 


± 


0.03 


0. 


.76 


± 





.02 


0.04 


5 


1 


-1.0862 


-1.0670 


Nd 


60 


1.46 


± 


0.03 


1. 


.45 


± 


0. 


.01 


0.05 


46 


4 


-0.3723 


-0.5163 


Sm 


62 


0.95 


± 


0.04 


1. 


.00 


± 





.01 


0.05 


36 


5 


-0.7595 


-0.7592 


Eu 


63 


0.52 


± 


0.04 





.52 


± 


0. 


.01 


0.04 


14 


6 


-1.0424 


-1.0376 


Gd 


64 


1.06 


± 


0.02 


1. 


.11 


± 


0. 


.01 


0.05 


20 


7 


-0.5591 


-0.5546 


Tb 


65 


0.31 


± 


0.03 





.28 


± 





.07 


0.10 


2 


8 


-1.2218 


-1.2526 


Dy 


66 


1.13 


± 


0.04 


1. 


.13 


± 


0. 


.02 


0.06 


13 


1 


-0.4437 


-0.4755 


Ho 


67 


0.49 


± 


0.02 





.51 


± 


0. 


.10 


0.10 


3 


9 


-1.0899 


-1.0862 


Er 


68 


0.95 


± 


0.03 


0. 


.96 


± 


0. 


03 


0.06 


8 


10 


-0.6798 


-0.6832 


Tm 


69 


0.11 


± 


0.06 





.14 


± 





.02 


0.04 


3 


1 


-1.5086 


-1.4841 


Yb 


70 


0.94 


± 


0.03 


0. 


.86 


± 





.10 


0.10 


1 


1 


-0.7889 


-0.7783 


Lu 


71 


0.09 


± 


0.06 


0. 


.12 


± 


0. 


.08 


0.08 


1 


1 


-1.5100 


-1.5317 


Hf 


72 


0.77 


± 


0.04 





.88 


± 





.02 


0.03 


4 


11 


-1.0974 


-1.1675 



a odders 2003 

b number of lines used for the photospheric abundance 
c reference for the photospheric abundance 

d Estimates of the r-process only abundances N r of solar-system RE elements, based 
on the differences between total meteoritic abundances N met and "empirical" and "stel- 
lar" estimates of the s-process only abundances N s ; see text for explanation of these 
estimates. These meteoritic abundances (normalized to log -ZV(Si) = 6) can be translated 
to photospheric ones (normalized to loge(H) = 12) through loge = logiV + 1.54 

References. — 1. This paper; 2. Lawler et al. (2001a); 3. Lawler et al. (2009); 4. 
Den Hartog et al. (2003); 5. Lawler et al. (2006); 6. Lawler et al. (2001c); 7. Den 
Hartog et al. (2006); 8. Lawler et al. (2001b); 9. Lawler et al. (2004); 10. Lawler et al. 
(2008b); 11. Lawler et al. (2007) 
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Table 4. Rare Earth Abundances for BD+17 3248, CS 22892-052, and CS 31082-001 











BD+17 3248 






CS 22892-052 








CS 31082-001 






El 


Z 


loge 




a 


# a 


Ref b 


loge 




(7 


# a 


Ref b 


lc 






<J 


# a 


Ref b 


Ba 


56 


+0.48 ± 0.05 





.11 


4 


1 


-0.01 ± 0.06 


0. 


.12 


4 


1 












1 


La 


57 


-0.42 


± 


0.01 





.05 


15 


2 


-0.84 ± 0.01 





.05 


15 


2 


-0.62 


± 0.01 





.04 


9 


2 


Ce 


58 


-0.11 


± 


0.01 





.05 


10 


4 


-0.46 ± 0.01 


0. 


.05 


32 


i 


-0.29 


± 0.01 





.03 


38 


4 


Pr 


59 


-0.71 


± 


0.02 





.06 


18 


1 


-0.96 ± 0.02 





.07 


15 


1 


-0.79 


± 0.01 


0. 


.07 


27 


1 


Nd 


60 


-0.09 


± 


0.01 





.06 


57 


5 


-0.37 ± 0.01 





.06 


37 


5 


-0.15 


± 0.01 


0. 


.08 


68 


1 


Sin 


(52 


-0.34 


± 


0.01 





.05 


72 


6 


-0.61 ± 0.01 


0. 


.07 


55 


6 


-0.42 


± 0.01 





.04 


67 


1 


Eu 


63 


-0.68 


± 


0.01 





.04 


9 


1 


-0.95 ± 0.01 





.02 


9 


1 


-0.72 


± 0.01 





.03 


7 


1 


Gd 


61 


-0.14 


± 


0.01 





.04 


11 


7 


-0.42 ± 0.01 





.07 


32 


7 


-0.21 


± 0.01 





.05 


32 


1 


Tb 


65 


-0.91 


± 


0.02 





.05 


5 


8 


-1.13 ± 0.01 


0. 


.04 


7 


9 


-1.01 


± 0.01 


0. 


.04 


9 


1 


Dy 


66 


-0.04 


± 


0.01 





.05 


28 


1 


-0.26 ± 0.01 


0. 


.06 


29 


1 


-0.07 


± 0.01 





.05 


35 


1 


Ho 


67 


-0.70 


± 


0.02 





.05 


11 


10 


-0.92 ± 0.01 





.02 


13 


10 


-0.80 


± 0.03 


0. 


.09 


12 


1 


Er 


68 


-0.25 


± 


0.01 





.04 


17 


11 


-0.48 ± 0.01 


0. 


.04 


21 


11 


-0.30 


± 0.01 





.04 


19 


11 


Tm 


69 


-1.12 


± 


0.02 





.05 


6 


1 


-1.39 ± 0.02 


0. 


.04 


6 


1 


-1.15 


± 0.02 





.06 


7 


1 


Yb 


70 


-0.27 


± 


0.10 





.10 


1 


1 


-0.55 ± 0.10 





.10 


1 


1 


-0.41 


± 0.10 


0. 


.10 


1 


1 


Lu 


71 














1 










1 












1 


Hf 


72 


-0.57 


± 


0.03 





.08 


6 


2 


-0.88 ± 0.01 


0. 


.04 


8 


2 


-0.72 


± 0.01 





.04 


10 


2 



a number of lines used for the stellar abundance 

b reference for the stellar abundance; these are cited at the end of Table [5] 
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Table 5. Rare Earth Abundances for HP 115444 and HP 221170 



HD 115444 HD 221170 



El 


Z 


loge 


a 


# a 


Ref b 


loge 


a 


#a 


Ref b 


Ba 


56 


-0.73 ± 0.04 


0.08 


4 


1 


+0.18 ± 0.05 


0.11 


4 


1 


La 


57 


-1.44 ± 0.02 


0.05 


8 


2 


-0.73 ± 0.01 


0.06 


36 


3 


Ce 


58 


-1.06 ± 0.01 


0.07 


26 


4 


-0.42 ± 0.01 


0.04 


37 


4 


Pr 


59 


-1.57 ± 0.02 


0.06 


10 


1 


-1.00 ± 0.02 


0.07 


21 


1 


Nd 


60 


-1.02 ± 0.01 


0.08 


37 


5 


-0.35 ± 0.01 


0.08 


63 


3 


Sm 


62 


-1.26 ± 0.01 


0.07 


67 


6 


-0.66 ± 0.01 


0.07 


28 


3 


Eu 


63 


-1.64 ± 0.02 


0.04 


8 


1 


-0.89 ± 0.03 


0.07 


7 


1 


Gd 


64 


-1.08 ± 0.01 


0.07 


29 


7 


-0.46 ± 0.04 


0.14 


11 


3 


Tb 


65 


-1.84 ± 0.04 


0.08 


3 


1 


-1.21 ± 0.03 


0.08 


8 


3 


Dy 


66 


-1.00 ± 0.01 


0.07 


24 


1 


-0.29 ± 0.01 


0.06 


25 


1 


Ho 


67 


-1.61 ± 0.01 


0.04 


9 


10 


-0.97 ± 0.02 


0.07 


8 


3 


Er 


68 


-1.22 ± 0.02 


0.07 


15 


11 


-0.47 ± 0.02 


0.08 


14 


11 


Tm 


69 


-2.06 ± 0.02 


0.04 


5 


1 


-1.39 ± 0.03 


0.06 


6 


1 


Yb 


70 


-1.43 ± 0.10 


0.10 


1 


1 


-0.48 ± 0.10 


0.10 


1 


1 


Lu 


71 








1 








1 


Hf 


72 


-1.51 ± 0.01 


0.03 


4 


2 


-0.84 ± 0.03 


0.11 


10 


2 



a numbcr of lines used for the stellar abundance 
b reference for the stellar abundance 

References. — 1. This paper; 2. Lawler et al. (2007); 3. Ivans et al. (2006); 4. Lawler 
ct al. (2009); 5. Den Hartog ct al. (2003); 6. Lawler et al. (2006); 7. Den Hartog et al. 
(2006); 8. Cowan ct al. (2002); 9. Sneden ct al. (2003); 10. Lawler et al. (2004); 11. 
Lawler et al. (2008b); 
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Table 6. Pr II Line Abundances 



A 


X 


log gf 


logs'/ !oge !oge 


loge 


loge 


log e log e 


A 


eV 


Li07 


IvOl BD+17 3248 


CS 22892-052 


CS 31082-001 


HD 115444 HD 221170 



3964.82 


0.055 


+0.069 


+0.121 




-0.75 


-1.00 






-1.03 


3965.26 


0.204 


+0.204 


+0.135 




-0.85 


-1.03 






-1.03 


4004.70 


0.216 


-0.250 








-0.90 


-0.71 






4015.39 


0.216 


-0.362 










-0.84 






4039.34 


0.204 


-0.336 










-0.76 




-0.98 


4044.81 


0.000 


-0.293 






-0.65 


-0.92 


-0.72 




-0.88 


4062.81 


0.422 


+0.334 






-0.59 


-0.83 


-0.63 




-0.83 


4096.82 


0.216 


-0.255 






-0.75 




-0.81 






4118.46 


0.055 


+0.175 








-0.85 


-0.68 






4141.22 


0.550 


+0.381 






-0.80 


-1.04 


-0.86 




-1.08 


4143.12 


0.371 


+0.604 


+0.609 




-0.68 




-0.71 


-1.49 




4164.16 


0.204 


+0.170 


+0.160 




-0.75 


-1.00 


-0.84 




-1.05 


4179.40 


0.204 


+0.459 


+0.477 




-0.58 


-0.98 


-0.79 


-1.49 


-0.88 


4189.48 


0.371 


+0.431 


+0.382 




0.72 


-1.02 


-0.86 


-1.64 


-1.03 


4222.95 


0.055 


+0.235 


+0.271 


+0.71 


-0.70 


-1.00 


-0.74 


-1.61 


-1.00 


4405.83 


0.550 


-0.062 


-0.037 








-0.71 






4408.82 


0.000 


+0.053 


+0.179 




-0.70 




-0.71 


-1.53 


-0.94 


4413.77 


0.216 


-0.563 










-0.73 






4429.13 


0.000 


-0.495 






-0.70 


-1.02 


-0.78 


-1.59 


-1.03 


4429.26 


0.371 


-0.048 


0.103 














4449.83 


0.204 


-0.261 


-0.174 




-0.70 




-0.76 




-0.97 


4496.33 


0.055 


-0.368 


-0.268 




-0.72 


-0.90 


-0.76 


-1.49 


-0.97 


4496.47 


0.216 


-0.762 
















4510.15 


0.422 


-0.007 


0.023 


+0.78 


-0.72 




-0.86 


-1.64 


-1.02 


5129.54 


0.648 


0.134 










-0.81 




-1.01 


5135.15 


0.949 




+0.008 








-0.91 




-1.03 


5173.91 


0.967 


+0.359 


+0.384 








-0.86 






5219.05 


0.795 


(+0.405) a 


-0.053 


+0.81 












5220.11 


0.795 




+0.298 




-0.72 


-1.00 


-0.89 


-1.59 


-1.08 


5259.73 


0.633 




+0.114 


+0.78 


-0.70 


-0.91 


-0.86 


-1.64 


-1.08 


5292.62 


0.648 


-0.269 


-0.257 








-0.83 




-1.05 


5322.77 


0.482 


-0.123 


-0.319 


+0.74 






-0.81 




-1.05 


5352.40 


0.482 


-0.739 

















107 note that this is a blended line in their spectrum; we used the log gf from IvOl. 
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Table 7. Dy II Line Abundances 



A 


X 




loge 


loge 


loge 


loge 


loge 


loge 


A 


eV 







BD+17 3248 


CS 22892-052 


CS 31082-001 


HD 115444 


HD 221170 


3407.80 


0.000 


+0.18 




-0.04 


-0.20 


-0.02 


-1.04 




3413.78 


0.103 


-0.52 




-0.01 


-0.23 


-0.07 


-1.01 




3434.37 


0.000 


-0.45 


+ 1.00 


-0.09 


-0.23 


-0.05 


-0.90 




3454.32 


0.103 


-0.14 


+ 1.20 


-0.09 


-0.33 


-0.17 




-0.28 


3456.56 


0.589 


-0.11 




-0.01 


-0.35 


-0.15 


-0.90 


-0.28 


3460.97 


0.000 


-0.07 




-0.19 


-0.28 


-0.13 


-1.04 


-0.28 


3523.98 


0.538 


+0.42 




-0.10 




-0.12 


-1.04 




3531.71 


0.000 


+0.77 


+ 1.20 


-0.02 


-0.23 


-0.03 


-1.11 


-0.18 


3534.96 


0.103 


-0.04 




-0.11 


-0.35 


-0.10 






3536.02 


0.538 


+0.53 


+ 1.10 


-0.07 


-0.35 


-0.05 


1.11 


-0.23 


3546.83 


0.103 


-0.55 


+ 1.23 


-0.01 


-0.23 


-0.07 


-0.91 


-0.30 


3550.22 


0.589 


+0.27 




-0.04 


-0.35 


-0.17 


-1.06 


-0.31 


3551.62 


0.589 


+0.02 






-0.32 


-0.15 






3563.15 


0.103 


-0.36 




-0.06 


-0.30 


-0.13 


-1.05 


-0.30 


3630.24 


0.538 


+0.04 




-0.05 


-0.30 


-0.09 


-0.94 


-0.27 


3630.48 


0.925 


-0.66 








-0.07 






3694.81 


0.103 


-0.11 


+ 1.11 


-0.06 


-0.28 


-0.08 


-1.03 


-0.25 


3747.82 


0.103 


-0.81 




-0.05 


-0.28 


-0.08 


-0.92 


-0.34 


3757.37 


0.103 


-0.17 




-0.05 


-0.31 


-0.08 


-1.01 


-0.28 


3788.44 


0.103 


-0.57 




-0.02 


-0.23 


-0.07 


-0.96 


-0.35 


3944.68 


0.000 


+0.11 




-0.06 


-0.20 


+0.00 


-1.06 


-0.30 


3978.56 


0.925 


+0.22 






-0.30 








3983.65 


0.538 


-0.31 


+ 1.08 


-0.02 


-0.28 


-0.07 


-0.97 


-0.26 


3996.69 


0.589 


-0.26 


+ 1.10 


-0.03 


-0.26 


-0.08 


-0.96 


-0.36 


4011.29 


0.925 


-0.73 


+ 1.14 






-0.08 






4014.70 


0.927 


-0.70 






-0.23 


-0.03 






4041.98 


0.927 


-0.90 








-0.05 






4050.57 


0.589 


-0.47 




-0.03 


-0.23 


-0.05 


-1.01 


-0.38 


4073.12 


0.538 


-0.32 


+ 1.10 


-0.02 


-0.23 


-0.07 


-1.06 


-0.41 


4077.97 


0.103 


-0.04 


+ 1.17 


-0.01 


-0.16 


-0.03 


-0.94 


-0.28 


4103.31 


0.103 


-0.38 


+ 1.17 


+0.11 


-0.15 


+0.02 


-0.91 


-0.13 


4124.63 


0.925 


-0.66 




-0.04 


-0.26 


-0.02 




-0.30 


4409.38 


0.000 


-1.24 




+0.03 




-0.02 




-0.28 


4449.70 


0.000 


-1.03 


+ 1.12 


+0.06 


-0.22 


-0.05 


-0.96 


-0.38 


4620.04 


0.103 


-1.93 








+0.05 




-0.28 


5169.69 


0.103 


-1.95 








-0.05 




-0.33 



Table 7— Continued 



A x log 5/ 


loge 


loge 


loge 


loge 


log e log e 


A eV 





BD+17 3248 


CS 22892-052 


CS 31082-001 


HD 115444 HD 221170 
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Table 8. Tm II Line Abundances 



A 


X 


logg/ 


loge 


loge 


loge 


loge 


loge 


loge 


A 


eV 







BD+17 3248 


CS 22892-052 


CS 31082-001 


HD 115444 


HD 221170 


3240.23 


0.029 


-0.80 








-1.09 






3462.20 


0.000 


+0.03 


+0.18 


-1.20 


-1.37 


1.24 


-2.13 


-1.37 


3700.26 


0.029 


-0.38 


+0.13 


-1.06 


-1.34 


-1.11 


-2.04 


-1.35 


3701.36 


0.000 


-0.54 


+0.10 


-1.10 


-1.42 


1.17 


-2.04 


-1.35 


3795.76 


0.029 


-0.23 




-1.13 


-1.44 


-1.19 


-2.07 


-1.52 


3848.02 


0.000 


-0.14 




-1.12 


-1.37 


-1.19 


-2.02 


-1.37 


3996.51 


0.000 


-1.20 




-1.10 


-1.37 


-1.07 




-1.40 
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Table 9. Abundance Sensitivities to Parameter Changes 



parameter= 
change 13 


+150 


logs 
+0.5 


v t 
+0.5 


[M/H] 
+0.5 


scat a 
yes 


model 
MARCS C 


La 


+0. 


.10 


+0. 


.16 


-0. 


.02 


+0. 


,03 


-0. 


.07 


-0 


.01 


Ce 


+0. 


09 


+0. 


14 


-0. 


.02 


+0. 


,04 


-0. 


.05 


-0 


.01 


Pr 


+0. 


.11 


+0. 


.14 


-0. 


.02 


+0. 


,04 


-0. 


.05 


0. 


,00 


Eu 


+0, 


.11 


+0. 


.16 


-0. 


.03 


+0. 


,03 


-0. 


.05 


+0 


.01 


Dy 


+0. 


.10 


+0. 


.12 


-0. 


.05 


+0, 


,03 


-0. 


.10 


-0. 


.01 


Er 


+0. 


.10 


+0. 


.11 


-0. 


.08 


+0, 


,03 


-0, 


.12 


-0. 


.01 


Tm 


+0. 


09 


+0, 


13 


-0. 


.03 


+0 


,04 


-0. 


.10 


0, 


,00 


Yb 


+0. 


.08 


+0. 


.05 


-0. 


,20 


-0, 


.05 


-0, 


23 


-0 


.06 


<> 


+0, 


.10 


+0. 


13 


-0. 


.06 


+0, 


.02 


-0, 


,10 


-0 


,01 


a 


0. 


01 


0. 


04 


0. 


,06 


0, 


.03 


0, 


06 


0. 


,03 



a Continuum source function computed with scattering. 

b Change of a model parameter from a baseline model taken from 
Kurucz (1998) grid computed with parameters T e ff = 4750 K, 
logs = 1-5, Vt = 2.0 km s _1 , [M/H] = —2.5, and no correction for 
scattering opacity in the continuum source function. 



c Gustafsson et al. (2008). 
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Table 10. Measured hfs A Constants for Pr II Levels of Interest in This Investigation 



Energy" 


Energy 


J 


hfs A 


Ref. 


cm -1 


cm -1 




0.001 cnr 1 




0.00 


0.000 


4 


-7.962 ± 0.013 


b 








-7.3 ± 0.9 


c 








-7.3 


d 


441.95 


442.079 


5 


63.721 ± 0.070 


b 








63.8 ± 1.1 


c 








63.9 


d 


1649.01 


1649.092 


6 


54.498 ± 0.090 


b 








53.8 ± 1.8 


c 








54.5 


d 


1743.72 


1743.776 


5 


-1.478 ± 0.030 


b 








-1.6 ± 0.5 


c 








-1.3 


d 



References. — a. Martin et al. (1978); b. Rivest 
et al. (2002); c. IvaOl; d. Ginibrc (1989); e. Li 
et al. (2000b); f. Li et al. (2000a); g. Ma et al. 
(1999) 

Note. — This table is presented in its entirety as 
a formatted ASCII file in the electronic edition of 
the Astrophysical Journal Supplement. A portion 
is shown here for guidance regarding its form and 
content. 
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Table 11. Hyperfine Structure Line Component Patterns for Pr II 



Wavenumber a \ a ir a F upp Fi ow Component Component Strength 

Position 13 Position 13 



cm 


-l 


A 










cm 1 


A 






25467. 


.565 


3925 


.4515 


6. 


.5 


6 


.5 


0.32402 


-0.049944 


0. 


.23932 


25467. 


.565 


3925 


.4515 


6. 


.5 


5 


.5 


0.27227 


-0.041967 


0. 


.01994 


25467. 


.565 


3925 


.4515 


5. 


.5 


6 


.5 


0.16516 


-0.025458 


0. 


.01994 


25467. 


.565 


3925 


.4515 


5 


.5 


5 


.5 


0.11341 


-0.017481 





,17164 


25467. 


.565 


3925 


.4515 


5. 


.5 


4 


.5 


0.06962 


-0.010731 





.03064 


25467. 


.565 


3925 


.4515 


4. 


.5 


5 


.5 


-0.02101 


0.003239 


0. 


.03064 


25467. 


.565 


3925 


.4515 


4, 


.5 


4 


.5 


-0.06480 


0.009989 





,12121 


25467. 


.565 


3925 


.4515 


4. 


.5 


3 


.5 


-0.10063 


0.015512 


0. 


.03333 


25467. 


.565 


3925 


.4515 


3. 


.5 


4 


.5 


-0.17479 


0.026941 


0. 


.03333 


25467. 


.565 


3925 


.4515 


3 


.5 


3 


.5 


-0.21062 


0.032464 





,08571 


25467. 


.565 


3925 


.4515 


3 


.5 


2 


.5 


-0.23848 


0.036760 


0. 


.02910 


25467. 


.565 


3925 


.4515 


2 


.5 


3 


.5 


-0.29616 


0.045650 


0. 


,02910 


25467. 


.565 


3925 


.4515 


2. 


.5 


2 


.5 


-0.32402 


0.049945 





.06349 


25467. 


.565 


3925 


.4515 


2. 


.5 


1 


.5 


-0.34393 


0.053014 





,01852 


25467. 


.565 


3925 


.4515 


1 


.5 


2 


.5 


-0.38512 


0.059364 


0. 


,01852 


25467. 


.565 


3925 


.4515 


1. 


.5 


1 


.5 


-0.40503 


0.062432 


0, 


.05556 



a Center-of-gravity value 

b Rclative to the center-of-gravity value 

c Normalizcd to 1 for the whole transition 

Note. — This table is presented in its entirety as a formatted ASCII file in the 
electronic edition of the Astrophysical Journal Supplement. A portion is shown here 
for guidance regarding its form and content. 
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Table 12. Isotopic and Hyperfine Structure Line Component Patterns for Yb II 



Wavenumber a A a i r a F upp Fi ow Component Component Strength 

Position b Position 13 Isotope 
cm -1 A cm -1 A 



30392 


.23 


3289. 


.367 


1 


.5 





.5 


0. 


.09884 


-0 


.010697 


0. 


.00130 


168 


30392 


.23 


3289. 


.367 


1 


.5 





.5 


0. 


.06745 


-0 


.007300 


0, 


.03040 


170 


30392 


.23 


3289, 


.367 


1 


.5 





.5 


0, 


.01878 


-0 


.002033 


0. 


.21830 


172 


30392 


.23 


3289. 


.367 


1 


.5 





.5 


-0, 


,01971 





.002134 


0. 


,31830 


174 


30392 


.23 


3289. 


.367 


1 


.5 





.5 


-0. 


,05657 





.006123 


0. 


,12760 


176 


30392 


.23 


3289. 


,367 


2 


.0 


1 


.0 


-0, 


.03402 





.003682 


0, 


,08925 


171 


30392 


.23 


3289. 


,367 


1 


.0 


1 


.0 


-0 


.09253 





.010015 


0. 


,01785 


171 


30392 


.23 


3289. 


,367 


1 


.0 





.0 


0. 


32919 


-0 


.035629 


0. 


.03570 


171 


30392 


.23 


3289. 


,367 


4 


.0 


3 


.0 


0. 


.12828 


-0 


.013884 


0. 


.06049 


173 


30392 


.23 


3289, 


,367 


3 


.0 


3 


.0 


0. 


12201 


-0 


.013206 


0. 


,02614 


173 


30392 


.23 


3289. 


367 


3 


.0 


2 


.0 


-0. 


.22795 





.024673 


0. 


,02091 


173 


30392 


.23 


3289. 


,367 


2 


.0 


3 


.0 


0. 


.16844 


-0 


.018231 


0. 


.00747 


173 


30392 


.23 


3289, 


,367 


2 


.0 


2 


.0 


-0. 


.18152 





.019647 


0. 


.02614 


173 


30392 


.23 


3289. 


,367 


1 


.0 


2 


.0 


-0. 


.12622 





.013661 


0. 


.02016 


173 


27061 


.82 


3694. 


,192 





.5 





.5 


0. 


10060 


-0 


.013732 


0. 


.00130 


168 


27061 


.82 


3694 


,192 





.5 





.5 


0. 


07469 


-0 


.010196 


0, 


.03040 


170 


27061 


.82 


3694. 


,192 





.5 





.5 


0. 


.02054 


-0 


.002804 


0. 


.21830 


172 


27061 


.82 


3694 


,192 





.5 





.5 


-0. 


,02200 





.003003 


0, 


,31830 


174 


27061 


.82 


3694. 


,192 





.5 





.5 


-0. 


.06261 





.008547 


0. 


,12760 


176 


27061 


.82 


3694. 


,192 


1 


.0 


1 


.0 


-0. 


.03284 





.004483 


0. 


,07140 


171 


27061 


.82 


3694. 


,192 


1 


.0 





.0 


0. 


.38888 


-0 


.053087 


0. 


.03570 


171 


27061 


.82 


3694. 


,192 





.0 


1 


.0 


-0, 


.10305 





.014067 


0. 


.03570 


171 


27061 


.82 


3694 


,192 


3 


.0 


3 


.0 


0. 


,12312 


-0 


.016807 


0. 


,04182 


173 


27061 


.82 


3694. 


,192 


3 


.0 


2 


.0 


-0. 


.22685 





.030968 


0. 


.05227 


173 


27061 


.82 


3694. 


.192 


2 


.0 


3 


.0 


0. 


,18128 


-0 


.024746 


0. 


.05227 


173 


27061 


.82 


3694. 


,192 


2 


.0 


2 


.0 


-0. 


.16869 





.023029 


0. 


,01494 


173 



a Ccnter-of-gravity value 

b Relative to the center-of-gravity value 

c Normalized to 1 for the whole transition 

Note. — This table is presented as a formatted ASCII file in the electronic edition of the 
Astrophysical Journal Supplement. 



